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The trend in electronic circuit design is toward increasing power dis- 
sipation density. The performance and reliability of increasing number of 
electronic systems are now seriously threatened by thermal effects, so that 
it is necessary to reappraise the relevant thermal design procedures. This 
paper examines natural convection cooling and concerns the prediction 
of maximum temperatures of electronic cabinets containing arrays of 
vertically oriented circuit cards with unequal power dissipation levels. By 
means of a vertical channel model, the effects of channel spacing, channel 
height, and power dissipation level are assessed with emphasis on asym- 
metric powering of the channel walls. Methods are indicated for rapid 
evaluation of maximum temperatures and optimum channel spacing with 
asymmetric heating. The present results show that asymmetry reduces the 
thermal performance of the channel. Consequently, the power dissipation 
on the channel walls should be made nearly equal. 

I. INTRODUCTION 

For the past several years, the trend in electronic equipment design 
has been toward ever-increasing circuit speeds. Now it is common for 
response times of telecommunication equipment to be specified in 
terms of nanoseconds while, in high-performance data processing sys- 
tems, the required response times are given in the picosecond range. 

The increase in circuit speed is facilitated by integration of circuit 
functions, circuit miniaturization, and higher-density packaging. Al- 
though the miniaturization of circuits results in decreased power dis- 
sipation per circuit, the power generation per unit volume, which is the 
important parameter in determining the circuit temperature, is actually 
increased due to the much higher packaging densities. The thermal 
problem is also compounded by the lower operating temperature re- 
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quirements of integrated circuits. This gives rise to a challenging under- 
taking in thermal design. And yet a poor thermal design could pos- 
sibly lead to complete failure or unreliable performance of equipment. 
Thus, it is imperative that a sound thermal design be initiated, and this 
at the earliest possible stage of system planning. 

The thermal design of a modern electronic system is based on a ra- 
tional selection of a cooling option followed by a thoughtful design con- 
sideration. A design should not only be practicable, but also economi- 
cal, serviceable, reliable, and compatible with other system compo- 
nents. It is the purpose of this paper to present design data on natural 
convection cooling of modern electronic systems. The emphasis is on 
the prediction of maximum temperatures in equipment. Although na- 
tural convection is the oldest method used in electronics cooling, the 
method is still employed extensively in new generations of equipment 
both in the communications and the data processing fields. In many 
cases, this method is used in conjunction with one or more of the newer 
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Fig. 1 — Schematic diagram of electronic equipment. 
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Fig. 2 — Two-dimensional vertical channel model. 



methods such as forced flow and the application of heat pipes and 
liquid cooling. 

The present paper concerns equipment wherein circuit boards con- 
taining heat-dissipative components are mounted vertically in channel- 
like fashion as shown in Fig. 1. Since the circuit boards are not iden- 
tically powered, the average heat fluxes from the columns are different, 
giving rise to asymmetrical air flow. In each channel, the heat transfer 
may be modeled by that in a vertical, two-dimensional, smooth-walled 
channel as shown in Fig. 2. The channel walls are treated as uniformly 
heated. However, the two wall heat fluxes need not be identical. 

A number of investigations have been reported in the literature 
concerning free convection in vertical channels. Ostrach 1 solved the 
combined free (natural) and forced convection problem in a vertical 
channel in fully developed flow with symmetrical uniform heat flux and 
internal heat generation. By fully developed flow is meant a situation 
in which the fluid velocity is invariant in the direction of flow. Engel 
and Mueller 2 investigated the effect of nonisothermal channel walls by 
assuming constant heat fluxes. Their results, however, are limited to 
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symmetrical wall heating. Lauber and Welch 3 considered the combined 
free and forced convection between vertical flat plates which are heated 
at uniform but different heat fluxes. Their work is confined to fully 
developed flow. 

Two standard references widely applied to heat transfer in electronic 
equipment are the experimental work of Elenbaas 4 and the numerical 
solution of Bodoia and Osterle. 5 Inasmuch as the above references deal 
with a channel whose walls are at a constant temperature, two factors 
are therefore ignored in their results: the effect of nonisothermal sur- 
faces and asymmetric heating. These effects are considered in the 
present paper. 

It is to be noted that, as in the references cited above, radiative trans- 
fer of heat is ignored in this paper. In some electronic systems, radiation 
may be a significant factor in heat removal. Inasmuch as temperatures 
of heat sources are of interest, the neglect of radiation usually results 
in conservative (too high) estimates of -the temperatures. 

II. ANALYSIS 

Referring to Fig. 2 which shows a two-dimensional, straight vertical 
channel, let the channel height be I and width be b. The channel walls 




g/3qb° 



Fig. 3 — Relation between dimensionless maximum wall temperature rise and the 
dimensionless parameter L (from Aung, et al. 6 ). 
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are heated causing fluid to rise between them by free convection. The 
thermal conditions of the walls are characterized as uniform heat flux 
(UHF), but the individual values for the two walls need not be the 
same. The fluid that enters the channel at the ambient temperature To 
is assumed to have a flat velocity profile Wo- The equations expressing 
conservation of mass, momentum, and energy for free convection in 
the constant-property (except for density in the buoyancy term) 
laminar flow may be written respectively as: 

»i+-S='0-7^ + « 5(r - r '>- (2) 

where p' is the pressure excess above the hydrostatic pressure. 
The boundary conditions are: 

For x = and < y < b : 

u = wo, T = T ; 

For y = and x ^ 0: 



u = 0, v = 0, k— = — q^; 



dy 



For y = b and x ^ 0: 



(4) 



u = 0, v = 0, k — = qi. 

At x = and x = £: p = po, where po is the hydrostatic pressure. 

The system of eqs. (1) to (4) may be cast in nondimensional form and 
then solved by numerical integration using a digital computer. Details 
are contained in Refs. 6 and 7. The latter references show that the nu- 
merical results agree well with experimental data. It may be pointed 
out that, for a channel whose height is large compared to the spacing so 
that the so-called fully developed flow condition exists, the governing 
eqs. (1) to (3) may be simplified and the solution is then obtainable in 
closed form. This is indicated in Ref. 8. In the present paper, these re- 
sults will be applied to indicate the effect of operating parameters on 
maximum wall temperatures. The cooling fluid considered is exclusively 
air. 
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III. RESULTS 

When the channel walls are individually heated at uniform heat 
fluxes, the quantity that is of design interest is generally the maximum 
temperature rise on the hotter wall. This quantity may be obtained 
from the general results shown in Fig. 3. It is clear from the latter that 
the maximum temperature increase depends on the parameter L 
which in turn depends on the quantity l/qb 5 . When the latter quantity 
is large, fully developed flow is approached 6 and the desired informa- 
tion on the maximum temperature increase is given by a rather simple 
expression which follows. 

3.1 Results at Large l/qb h 

At large values of the quantity t/qb 5 , the flow in the channel exhibits 
interesting characteristics whereby the velocity distribution across the 
channel remains unaltered with axial distance. All axial temperature 
variations are also given by linear relations. Typical velocity and tem- 
perature distributions in this situation are given in Fig. 4. This is the 
so-called fully developed flow situation. Clearly, it prevails most com- 
monly in a channel whose height I is large compared to its width. In 
a channel with developing flow, typical velocity and temperature dis- 
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Fig. 4 — Typical temperature and velocity distributions in a channel having nearly 
fully developed flow. 
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Fig. 5 — Typical temperature and velocity distributions in a channel having de- 
veloping flow. 

tributions are as indicated in Fig. 5. In both Fig. 4 and Fig. 5, the quan- 
tity M is a dimensionless mass flow parameter which is related to L 
(see Refs. 6 and 8). 

Figure 3 shows that, as fully developed flow is approached, all curves 
coalesce into one. Based on the derivations given in Ref. 8, it can be 
shown that, in fully developed uniform heat flux channels with asym- 
metric heating, the maximum temperatures on the hotter and cooler 
walls become practically identical and may be evaluated using the 
equation : 

flmax.l = 0max.2 = 6.928o(Pr)" 1 / 2 (L) 1 ' 2 , 

L^b, all r H 



(5) 



For air at 120°F, the Prandtl number (Pr) is about 0.7. Evaluating the 
thermal physical properties for air in eq. (5) at a temperature of 120°F, 
we may rewrite eq. (5) as: 



ci - t = 7Va*. 2 - t = o.ieiofor-w-W)- 1 - 6 ' 

all r H 



for ^> 1.19X10", 



(6) 



In eq. (6), I is expressed in ft, q in W/ft 2 , and b in ft. Note that, if T max ,i 
alone is to be evaluated, eq. (G) is applicable at l/qb h ^ 4.76 X 10 7 . 



914 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1973 

The reason is apparent from Fig. 3 where at L ^ 0.2 the deviations 
among curves of d maxA at all values of r H may be neglected. 

3.2 Results at Small C/qb 5 

The counter-case to the situation just discussed is the one in which 
the channel height I is relatively short compared to its width b. Again, 
in this situation the maximum temperature rise on the channel walls 
may be calculated by a simple relation. Note that, to use this relation, 
which will be given below, it is not necessary that I be smaller than b. 
It is only necessary to fulfill the validity limit attached to the equation. 

For relatively short channels, the results in Ref. 6 indicate that the 
maximum temperature in each wall is independent of the other wall, 
and may be calculated using the result for a single vertical flat plate. 9 
The latter result is valid for the asymmetrically heated channel when 
L ^ 10" 3 . Thus, we may write: 

max = 2.05(L) 1 ' 5 at I ^ 10" 3 . (7) 

Again, inserting thermal physical properties of air at 120°F, the above 
may be written: 

T max -T = 8.66^/y ' s j 

at ^- 6 ^ 2.35 X 10 5 , all r^ | 

In eq. (8), (Tmax - T ) is in degrees Fahrenheit, I in feet, and q in 
watts per square foot of surface area. Equation (8), which gives the 
maximum temperature on either the hotter or cooler wall when q is 
appropriately replaced by qi or q 2l shows that the maximum tempera- 
ture on the hotter wall, subject to the attached condition, is indepen- 
dent of the channel spacing which is to be expected. It may also be 
noted for reference that 

{ ( 1±IE \ i (9) 

qi b* ~ V 2 ) qb' K J 

3.3 Results at Intermediate l/qb* 

If the proposed design is such that neither eq. (6) nor eq. (8) may be 
applied, then the maximum temperature increase in the channel can 
be evaluated using Fig. 3. Parametric curves can be generated from 
Fig. 3 to display the effect of various operating variables on the maxi- 
mum temperature rise on the hotter wall. Since the value of L depends 
on the channel spacing, the average heat flux, and the channel height, 
any one of these may be varied while others remain constant at typical 
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Fig. 7 — Effect of channel height on the maximum wall temperature rise. 
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Fig. 8 — Effect of average heat flux on the maximum wall temperature rise, 
operating values. The following conditions are chosen as an example: 



b = 0.0365 foot (0.4375 inch) 

I = 6 feet 

q = 5.75 watts per square foot. 



(10) 



The resulting curves are shown in Figs. 6 through 8 which can be used 
by the designer. To obtain Fig. 6, I and q are obtained from eqs. (10) 
and, corresponding to each pair of assigned values of b and r H , L is 
calculated. A value of mB x.i is then obtained from Fig. 3. The maximum 
temperature rise can then be obtained. By varying r H and b, the solid 
curves in Fig. 6 result. The range of validity of eq. (8) is also indicated. 
Figures 7 and 8 are obtained in the same manner. The dashed lines in 
Figs. 6 and 8 pertain to optimum spacing which will be discussed later. 
From Figs. 6 through 8, it can be seen that at any fixed maximum tem- 
perature rise the effect of asymmetry is to increase the necessary chan- 
nel spacing, decrease the channel height, or decrease the heat flux when 
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any one of these quantities is taken as the sole variable in the 
equipment. 

The effect of asymmetric wall heating may be more conveniently 
examined by replotting Fig. 3 in the manner of Fig. 9. In the latter, 
the maximum temperature increase on each wall of an asymmetric 
channel is normalized by the maximum temperature increase on a 
similar channel wherein the same average heat flux is distributed evenly 
on the two walls, so that r H = 1. Clearly, as the degree of asymmetry 
is increased so that r H is decreased without affecting the average heat 
flux, the result is that the maximum temperature on the hotter wall is 
raised while that on the cooler wall is depressed. Therefore, in equip- 
ment design, it is desirable to obtain equal heating on the two channel 
walls. This is also desirable from the standpoint of optimum spacing, 
as will be seen. 

With the help of Fig. 9, the maximum temperature increases in an 
asymmetric channel may be evaluated once the symmetric heating 
value is known. The latter may be obtained rapidly with the aid of the 
nomogram shown in Fig. 10 once I, b, and q are known, r H being 1 in 
Fig. 10. The use of the latter has been described in Ref. 7 but is repeated 
here for ease of reference. To use Fig. 10, first locate Points 1, 2, and 4. 
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Draw a line connecting Points 1 and 2. Extend this line to intersect 
with the Pivot Axis /3 to give Point 3, and with Pivot Axis a to give 
Point 6. Draw a line joining Points 3 and 4. The intersection of this 
line with Pivot Axis 7 gives Point 5. Draw a line connecting Points 5 
and 6. Read the a value at Point 6 which for the case illustrated is 2. 
The intersection of line 5-6 with the appropriate a-branch (a = 2 in 
this case) in the center of the nomogram (Point 7) gives the maximum 
card temperature rise over room ambient, which in this example is 
38°C or 68°F. It should be noted that each a-branch is made up of 
three segments: a left-side vertical line which corresponds to flow 
around a single flat plate, a right-side vertical line which corresponds to 
fully developed flow, and an inclined straight line joining the two verti- 
cal ones. Moreover, the asymptotes of the vertical lines coincide for all 
values of a. Hence, if the line joining Points 5 and 6 does not intersect the 
appropriate a-branch on its inclined segment, either the left or the 
right vertical segment should be chosen depending on whether the 5-6 
line passes above the inclined segment or below it. 

3.4 Optimum Spacing Including the Effect of Asymmetry 

Consideration may now be given to the overall performance of the 
channel as a heat-removing device. In electronic equipment the maxi- 
mum device junction temperature is usually specified. Since the junc- 
tion temperature is related to the wall temperature, the maximum wall 
temperature is therefore implied. In modern equipment, it is desirable 
to increase the packing density and hence the total power dissipation of 
the entire equipment without increasing the maximum wall tempera- 
ture beyond the allowable limit. As shown in the appendix, the task 
here is to select b so as to maximize the heat transfer. The latter value 
will be designated g op t and the corresponding spacing is called 6 op t. 
Design curves for finding 6 0P t and g 0P t will be given below. The emphasis 
here is on asymmetrically heated channels. A more detailed discussion 
of optimum spacing in a symmetric channel may be found in Ref. 7. 
Note that maximum power dissipation is realized only when the equip- 
ment is strictly operating at 6 0P t and g opt . If a different (smaller) spac- 
ing is used, then a different (smaller) heat flux must be employed to 
yield the same maximum temperature increase. For this purpose, de- 
sign curves such as those in Figs. 6 through 8 can be consulted. 

Following Bodoia, 10 it can be shown that the power dissipation is 
maximized if the channel parameters are so selected that the channel is 
operated at the point where the slope of a log-log plot of Nu versus Ra 
is one-half, where Nu is the Nusselt number and Ra the Rayleigh 
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number denned as: 

Nu = 



qb 



1 



Ra 



(Tmax.l - T%)k m «,i' 

= felZ r Q )6* = o.7 X •_. X I. 



From Fig. 3, a relation between Nu and Ra can be obtained. This is 
shown in Fig. 11. Calling the values of Nu and Ra for maximum power 
dissipation Nu 0P t and Ra 0P t, respectively, Table I may be constructed 
from Fig. 11. In Table I, E is an efficiency of heat transfer defined as 



E = 



V 



where 



ill) 



p = Nu pt/(Raopt) 1/2 

&opt 

In the above, h opt and 6 opt are optimum values of A ma x.i and b which give 
maximum power dissipation. 




Fig. 11 — Relation between Nusselt and Rayleigh numbers. 
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Table I — Optimum Rayleigh and Nusselt Numbers 



rii 


Raopt 


Nuopt 


E 





42 


0.43 


65% 


0.1 


51 


0.51 


70% 


0.5 


70 


0.73 


86% 


1.0 


135 


1.18 


100% 



From Table I, it is seen that the efficiency E varies in the same 
manner as r«. Asymmetry therefore decreases the efficiency of heat 
transfer. From the thermal standpoint all equipment should therefore 
be designed to yield symmetrical heating as closely as possible. 

Using values of Ra 0P t from Table I, it is possible to obtain optimum 
spacings for different maximum temperatures once the channel height 
is specified. In like manner, the necessary heat flux at optimum spacing 
corresponding to different prescribed maximum temperatures can also 
be obtained from Nu op t- Results have been obtained for a 6-foot chan- 
nel with r H = 0.0 and th = 0.5 and compared with the case r H = 1 
in Figs. 6 and 8 (indicated by the dashed line). It may be ascertained 
from these figures that the effect of asymmetric heating is to decrease 
6 oP t and q ov t- The net effect of asymmetry is to decrease the total power 
dissipation in the cabinet (see appendix). 
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Fig. 12 — Optimum spacing and optimum heat flux as a function of the heat flux 
ratio th. 
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Fig. 13.— Nomogram for evaluation of optimum spacing, r H = 1 (from Aung, 
Kessler, and Beitin 7 ). 



Using the values of Ra op t and Nu op t given in Table I, it is possible to 
obtain a general relation between optimum spacing 6 0P t or optimum 
heat flux g op t and r H , since 



>o P t 



(&opt) rn=l 



(Ra opt ) 1/4 
(Ra 0P t)# 4 -i 



and 



3 



opt 



Nu 



opt 



(&pt) 



fff-1 



(Nuopt) 



r H =l 



The results are indicated in Fig. 12. With the aid of this figure and the 
nomogram given in Fig. 13 which yields (6 op t)r H =i, the optimum spac- 
ing at specified r H) I, and maximum hotter wall temperature rise 
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(T max ,i — To) can be rapidly evaluated. Note that the maximum tem- 
perature rise to be entered in Fig. 13 is also the specified value of 

(T max ,i — To). 

IV. CONCLUDING REMARKS 

In the present paper, consideration has been given to the free or 
natural convection cooling of electronic cabinets containing arrays of 
vertical circuit cards with asymmetric heating. It was indicated that 
fully developed flow prevailed in a channel whose height is large com- 
pared to its width. In this case, the maximum temperature in the chan- 
nel can be calculated using the simple relation presented here. In con- 
trast, when the channel spacing is relatively large, the maximum tem- 
perature on each wall is independent of the conditions on the other 
wall but can again be calculated by a simple equation. 

In situations where neither of the above two approaches can be ap- 
plied, the maximum temperature in the channel can be evaluated 
rapidly using the graphical procedure outlined in this paper. The opti- 
mum channel spacing for maximizing the power dissipation is discussed 
and it is emphasized that to reduce adverse thermal effects channel 
walls should be powered as equally as possible. 

APPENDIX 

The total power dissipation in the cabinet per unit depth can be 
written as 

Power = 2tNq 

= 2mh(T taKXtl - To), 

where AT is the number of channels (see Fig. 1). If the thickness of the 
channel wall is neglected, we have 



Hence, 



»-f 



Power = 2tW(T ma , 1 - T ) j- 

h 
a 6 



for fixed (T maXil — T ), W, and I. Clearly, to obtain the largest pos- 
sible power one needs to obtain the maximum h (that is q) and the 
minimum b, the combination of which gives the specified (T maXil — T ). 
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NOTATIONS 

b clear channel or card spacing in feet (ft) 

c p specific heat of air in W-s/lbm-F 

g acceleration due to gravity in ft/s 2 

h Q/(T max ,i - To), average heat transfer coefficient, in W/ft 2 -F 

k thermal conductivity of air in W/ft-F 

Li (& 2 k)/(g/3gi& 5 ), dimensionless 

L (& 2 k)/(g/3#& 6 ), dimensionless 

I channel height in ft 

Pr Prandtl number denned by MC P /k, dimensionless 

p pressure in lbf/ft 2 

p hydrostatic pressure in lbf/ft 2 

p' (p - p ) in lb,/ft 2 

q heat flux on channel wall per unit surface area, in W/ft 2 

q (<?i + 92)/2, W/ft 2 

th 92/91, dimensionless 

T temperature in degrees Fahrenheit (F) 

To ambient temperature in degrees Fahrenheit 

u, v axial and transverse velocity in ft/s 

W width of cabinet 

x, y axial and transverse coordinates in ft, see Fig. 2 

jfl thermal expansion coefficient of air in 1/F 

n dynamic viscosity in lb m /ft-s 

v kinematic viscosity of air in ft 2 /s 

p density of air in lb m /ft 3 

dimensionless temperature rise denned in Fig. 3 

Subscripts 

1 refers to hotter wall 

2 refers to cooler wall 
max maximum value 

max, 1 maximum value on hotter wall 
max, 2 maximum value on cooler wall 
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